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We study the stability of static, spherically symmetric solutions to the Einstein equations with a 
Q \ scalar field as the source. We describe a general methodology of studying small radial perturbations 

of scalar-vacuum configurations with arbitrary potentials V(4>) , and in particular space-times with 
^vq throats (including wormholes), which are possible if the scalar is phantom. At such a throat, the 

effective potential for perturbations V c $ has a positive pole (a potential wall) that prevents a 
q ' complete perturbation analysis. We show that, generically, (i) V^ff has precisely the form required 

for regularization by the known S-deformation method, and (ii) a solution with the regularized 
potential leads to regular scalar field and metric perturbations of the initial configuration. The well- 
^j" known conformal mappings make these results also applicable to scalar-tensor and f(R) theories 

of gravity. As a particular example, we prove the instability of all static solutions with both normal 
and phantom scalars and V{4>) = under spherical perturbations. We thus confirm the previous 
results on the unstable nature of anti-Fisher wormholes and Fisher's singular solution and prove 
the instability of other branches of these solutions including the anti-Fisher "cold black holes" . 



1. Introduction 



Phantom field configurations have become a subject of particular interest since the discovery of the 
accelerated expansion of our Universe and its explanation in the framework of general relativity by 
the existence of dark energy (DE), a source of gravity of unknown nature with the pressure to density 
ratio w < — 1/3. Numerous observations lead to estimates of w around —1, which corresponds to a 
cosmo logical constant, but values smaller than —1 are still admissible and, moreover, preferable for 
describing an increasing acceleration. The most recent estimates read w = —1.10 ± 0.14 (lcr) [1] 
(according to the 7-year WMAP data) and w = — 1.0691o;q92 [2] (mainly from data on type la 
supernovae from the SNLS3 sample). 

The possible existence of phantom DE is certainly only one of potentially viable explanations of 
the observations (numerous DE models are described, e.g., in the reviews [3-7]), but if we take it 
as a working hypothesis, it is natural to expect that there are manifestations of DE in local objects 
and phenomena. The simplest of them can be described by static, spherically symmetric solutions 
to Einstein-scalar equations where the scalar field has an unusual sign of kinetic energy (a phantom 
scalar, by definition). In the case of a massless scalar it is the so-called anti-Fisher solution, a 
phantom analogue of Fisher' solution [8] for an ordinary minimally coupled massless scalar field. 
The anti-Fisher solution consists of three branches, one of them, termed Branch C in this paper, 
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represents wormholes [9,10], the others, A and B, also have properties of interest, and some of them 
describe so-called cold black holes [11,12], which possess horizons of infinite area and zero Hawking 
temperature; in all such solutions there are throats (minima of the spherical radius), but, unlike 
wormholes, here beyond the throat one does not find an asymptotically flat region. 

A number of solutions for phantom scalar fields with nonzero self-interaction potentials are also 
known, among them are wormholes with flat and AdS asymptotic behaviours and a specific class 
of regular black holes called black universes [13,14]. These are configurations where a possible 
explorer, after crossing the event horizon, gets into an expanding universe instead of a singularity, 
and moreover, this expanding universe, being initially highly anisotropic, eventually isotropizes and 
approaches a de Sitter mode of expansion at late times. The stability of all such models under small 
perturbations is an important test of their possible viability [15]. 

The purpose of this paper is to present a methodology of studying small radial perturbations of 
scalar-vacuum configurations with any potential V(<f>) (see (1)), and in particular space-times with 
throats. The difficulty with the latter consists in the fact that the effective potentials V e s(x) for 
perturbations (not to be confused with the self-interaction potential V((f>)) always, in any gauge, 
possess a singularity at the throat, which prevents a complete perturbation analysis. It is for this 
reason that some previous stability studies of anti-Fisher solutions [11, 16] did not find any unstable 
mode whereas a numerical perturbation analysis of Shinkai and Hayward [17] revealed an instability 
of the simplest representative of this family of solutions, the Ellis massless wormhole [9, 10]. 

Gonzalez et al. [18], analyzing the stability of anti-Fisher wormholes, made a proper substitution 
in the perturbation equation (a special case of the so-called S-deformation [19,20]) and regularized 
the effective potential V c q. As a result, they found an exponentially growing mode with nonzero 
perturbation of the throat radius, thus showing that the anti-Fisher wormholes are unstable. We 
show here that a similar methodology can be applied to more general self-gravitating scalar field 
configurations including those with arbitrary self-interaction potentials V((f>) . To this end, we prove 
that, generically, (i) the effective potential V e s has precisely the form required for regularization by 
S-deformation and (ii) any solution of the transformed wave equation with a regularized potential 
leads to a regular perturbation of the background static configuration. 

As a particular example, we study the stability of all anti-Fisher solutions under spherically 
symmetric perturbations. We prove the instability of Branch A and B solutions and confirm the 
conclusions of [18] for Branch C solutions (wormholes). 

2. Perturbation equations 
2.1. Preliminaries 

Consider a self-gravitating, minimally coupled scalar field with an arbitrary self-interaction potential 
in general relativity. The Lagrangian is (up to a constant factor) 



where e = 1 for a normal scalar field with positive kinetic energy and e = — 1 for a phantom scalar 
field. Other notations are usual, the gravitational constant is absorbed in the definitions of <fi and 
V((f>). The field equations are 




(1) 



en 4> + v^ = o, 

where = dV/dcj). 



(2) 
(3) 
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The general spherically symmetric metric may be written in the form 

ds 2 = g^dx»dx u = e 2 ~<dt 2 - e 2a du 2 - e^dtt 2 , (4) 

where a, (3, 7 are functions of the radial coordinate u and the time coordinate t and dfl 2 = 
d9 2 + sin 2 9 dip 2 . We will also use the notation r = e 13 for the areal radius, such that 4irr 2 is 
the area of coordinate spheres u = const, t = const. There remains a coordinate freedom in the 
(u, t) subspace, which in general corresponds to choosing a reference frame preserving the spherical 
symmetry. In the static case, there is a reference frame such that there is no t-dependence, and 
then the coordinate freedom concerns the choice of the u coordinate. 

Curvature singularities for static, spherically symmetric space-times are entirely determined 
using the Kretschamnn scalar, which, for the metric (4), may be written as a sum of squares, 

K = 4K 2 + 8K% + SKl + 4Kl, (5) 

with 



la (3' 2 , (6) 

where the prime denotes d/du. The structure of Eq. (5) indicates that an infinite value of any Ki 
at some value of u implies the presence of a singularity at this u. All these expressions for Ki are 
invariant with respect to the choice of the u coordinate. 

We will assume that a certain static, spherically symmetric solution to Eqs. (2) and (3) is known 
and study its stability under small spherically symmetric perturbations. We thus consider, instead 
of <f>(u) , a perturbed unknown function 

<(>(u, t) = 4>{u) + 5cj)(u, i) 

and similarly for the metric functions a, (3, 7, where (f>(u) , etc., are taken from the static solutions. 

Preserving only linear terms with respect to time derivatives, we can write all the nonzero 
component of the Ricci tensor and the time-time component of the Einstein tensor as 

i?0 = e -2r (fi + 2/ 3) _ e -2 Q[ y, + y ( y _ Q , + ^ % (?) 

R\ = e ^a - e- 2a [ 7 " + 2/3" + 7 ' 2 + 2f3' 2 - a\i + 2/3% (8) 

R\ = + e- 2 ^(3 - e~ 2a [(3" + /3'( 7 ' - a' + 2(3')}, (9) 

Roi = 2[P' + Pf3' -a[3' -h% (10) 

G° Q = e- 2a [2/3" + /3'(3/3' - 2a')} - e" 2/3 , (11) 

where dots and primes denote d/dt and d/du, respectively. 
2.2. General form of the field equations 



The zero-order (i.e., static) scalar, Q , Q) , ( 2 ) components of Eqs. (3) and the Einstein equation 
0" + 0'(y + 2/3' - c/) =ee 2 % (12) 



G° = . . . are 
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7 " + 7 '( 7 ' + 2/3'-a') = -Ve 2a , (13) 

7" + 2(3" + 7 /2 + 2(3' 2 - a\i + 2/3') = -e0' 2 - V e 2a ; (14) 

- e 2a ~ 2/3 + (3" + /3'( 7 ' + 2/3' - a') = -V e 2a , (15) 

_ e 2«-2/3 + 2/ 3" + £'( 3/3 ' _ 2a ') = -I e 0' 2 - V" e 2a . (16) 

The first-order perturbed equations (scalar, R 01 — . . . , R\ = . . . , and Gq — . . . ) read 

e 2«-2i 5 _ 5( p" _ + 2(3' - a') - (j>'{5i + 25(3' - 5a') + e5( e 2a V<p) = 0, (17) 

6$' + (3' 5(3 - (3'5a - 7 '<5/3 = -\&l>'5<j>, (18) 

5( e 2a ~ 213 ) + e 2a - 27 5/3 - 5(3" - 5(3\i + 2(3' - a') 

- p'(5i + 25(3' - 5a') = 5{ e 2a V), (19) 

- 5(e 2Q - 2/3 ) + 5(3" + 6W - 2(3 '5a' - 2a' 5 P' = -e<j>' '5(f) - 5(Ve 2a ). (20) 

Eq. (18) may be integrated in t; since we are interested in time-dependent perturbations, we omit 
the appearing arbitrary function of u describing static perturbations and obtain 

8p' + 5/3(13' -7')- p'5a = -\e<i>'5<i>. (21) 



Let us note that we have two independent forms of arbitrariness: one is the freedom of choosing 
a radial coordinate u, the other is a perturbation gauge, or, in other words, a reference frame in the 
perturbed space-time, which can be expressed in imposing a certain relation for 5a, 5(3, etc. In 
what follows we will employ both kinds of freedom. All the above equations have been written in 
the most universal form, without coordinate or gauge fixing. 



2.3. Gauge 5/3 = 

This is technically the simplest gauge, in particular, it is convenient for considering usual black hole 
perturbations, but causes certain difficulties when applied to wormholes and other configurations 
with throats. The reason is that the assumption 5(3 = leaves invariable the throat radius, while 
perturbation must in general admit its time dependence [11]. This problem will be discussed below. 
With 5(3 = 0, Eq. (21) expresses 5a in terms of 5<p: 

2(3' 5a = e<t>'5<t). (22) 
Eq. (19) expresses 5^' — 5a' in terms of 5a and 5(f): 

$'{5i - 5a') = 2 e 2a - 2 ?5a - 5( e 2a V). (23) 
Substituting all this into (17), we obtain the wave equation 

e 2 «-27 5 _ S(j) " _ 5 0'( 7 ' + 2(3' - a') + U5(f) = 0, 

U - e 2 "{e(V - e- 2 ^ + + eV H ). (24) 

Before proceeding with a study of the wave equation, let us make sure that all the remain- 
ing Einstein equations hold as a consequence of (22), (17) and (19) and do not lead to any new 
restrictions. Consider the component (20) (the constraint equation). It now takes the form 

25a e 2a ~ w + 2(3' 5a' = e<f)'5<f) + 5(V e 2a ). (25) 
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This equation holds automatically owing to the zero-order equations and (22). Indeed, a substitution 
of 8a from (22) brings (25) to the form 

W - (ftp" + e 2a ~ 2/ y - e 2a V(f)' - e e 2c %/3'l = 0, (26) 

Now, substituting <fi" and j3" from (12) and (15), respectively, we see that all terms cancel, i.e., 
the equation does hold. Furthermore, the Einstein equation G\ — . . . holds as a consequence of 
(19) and (20); lastly, the equation G\ — ... holds due to the Bianchi identity V V G\ = and the 
corresponding property of the stress-energy tensor of the scalar field. 

So we can return to Eq. (24). Passing on to the "tortoise" coordinate x introduced according to 

du/dx = e 7 " (27) 

and changing the unknown function 8(f) (->■ if; according to 

8(f) = ip{x,t)e-P, <=> if>(x,t) =r8(f), (28) 

we reduce the wave equation to its canonical form, also called the master equation for radial per- 
turbations: 

$-if>xx + V eS (x)xjj = 0, (29) 

(the index x denotes d/dx), with the effective potential 

V eS (x) = e 2 ^ 2a [U + (3" + f3' 2 + p{1 - a')}. (30) 

This effective potential was previously obtained in other notations for e = —1 in [21]. A further 
substitution 

if){x, t) = y(x) e lujt , oo = const, (31) 

which is possible because the background is static, leads to the Schrodinger-like equation 

Vxx + - V cS (x)]y = 0. (32) 

If there is a nontrivial solution to (32) with Imu < satisfying some physically reasonable con- 
ditions at the ends of the range of u (in particular, the absence of ingoing waves), then the static 
system is unstable since 8(f) can exponentially grow with t. Otherwise our static system is stable 
in the linear approximation. Thus, as usual in such studies, the stability problem is reduced to a 
boundary-value problem for Eq. (32) — see, e.g., [11,18,22-26]. 

Note that all the above relations are written without fixing the background radial coordinate u . 

2.4. Gauge-invariant perturbations 

To be sure that we are dealing with real perturbations of the static background rather than purely 
coordinate effects, it is necessary to construct gauge-invariant quantities. 

Small coordinate transformations x a (-> x a + £ a in the (t, u) subspace can be written as 

t = t + At(t,u), u = u + Au(t,u), (33) 

where At and Au are supposed to be small. Any scalar quantity with respect to such transforma- 
tions, such as, e.g., (f>(t,u) acquires an increment: 

A(f) = (pAt + (f)'Au w (f)'Au (34) 
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in the linear approximation since both <fi and At are small. The quantity r, also being a scalar in 
the (t,u) subspace (a 2-scalar, for short), behaves in the same way. If there are perturbations 5(f) 
and Sr, the transformation (33) changes them as follows: 

5(f) h-> Jcj) = 5(f) + (j>'Au, 

5r (->■ 5r — 5r + r'Au. (35) 
It then follows that the combination 

V>i = r'5(j) - <f)'5r, (36) 

is invariant under the transformation (33), or gauge-invariant. Recall that the prime here denotes 
d/du in the background static configuration. 

One can notice that combinations constructed like (36) from any 2-scalars (for example, such as 
e^ and /3 = lnr , or two different linear combinations of <f) and r) are also gauge-invariant. Moreover, 
gauge-invariant is if)\ multiplied by any 2-scalar or any combination of background quantities which 
are known and fixed functions of u. 

The physical properties of perturbations must not depend on which gauge-invariant quantity if) 
is chosen to describe them. Meanwhile, with different if), the effective potentials will, in general, 
also be different. However, given a specific background configuration, in order that the theory be 
consistent, these different potentials should lead to the same perturbation spectrum. 

Due to gauge invariance of ip\ , equations that govern it may be written in any admissible gauge, 
in particular, the gauge 5(5 = , and Eq. (29) for if) then may be considered as a result of substituting 
if) = (r/r')tpi = r5(f) in a manifestly correct equation for ipi = r'5(f). 

On the similar problem regarding cosmological perturbations and the definition of the corre- 
sponding gauge invariants, see, e.g., the review [27]. 



2.5. Regularized potential near a throat 

The gauge 5/3 = (the same as 5r — 0) is suitable for describing the perturbations at any points 
except those where r' — 0: these are throats and other critical points of r{u). Indeed [11], putting 
5r = 0, we forbid perturbations of the throat radius, while there is no physical reason for that. 
Technically, this restriction manifests itself in a generically infinite value of the potential U(u) in 
Eq. (24) and consequently in V e s involved in the wave equation (29). Throats are only possible in 
the case e = — 1, and, provided U/r 2 < 1 at such a throat (u = u t h), the potential has there a 
wall of infinite height, with the generic behavior V^r ~ l/(u — u t h) 2 near the throat since we have 
there generically r'(u) ~ u — w t h- As a result, perturbations are actually independent at different 
sides of the throat, necessarily turn to zero at the throat itself, and we thus partly lose information 
on their possible properties. Such an incomplete treatment has led to a conclusion that anti-Fisher 
wormholes [9,10] (that is, wormhole solutions to Eqs. (2), (3) with V = and e = —1) were stable 
under spherically symmetric perturbations. A similar conclusion was made in [11] concerning such 
wormholes and cold black holes by using another (harmonic) gauge, 5a = 25(3 + 5^, which does not 
lead to a pole in the effective potential, but, as follows from our further consideration in this paper, 
this analysis was also incomplete. 

It could seem that the above difficulty only concerns the gauge 5/3 = 0. However, due to the 
gauge-invariant nature of Eq. (29) (to be verified below), it is clear that the problem is inherent to 
the background geometry itself, and the pole in the effective potential always emerges at a throat, 
if any. 
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A way of avoiding the restriction 5(3 («th) = is connected with the so-called S-deformations of 
the potential V e s. This method was used in [19,20] for transforming a partly negative potential 
to a positive-definite one in master equations for perturbations of higher-dimensional black holes. 
Using this method, Gonzalez et al. [18] transformed a singular potential to a nonsingular one for 
perturbations of the anti-Fisher wormholes and discovered the existence of an exponentially growing 
mode, showing that such wormholes are unstable. We will try to formulate a similar scheme suitable 
for the more general field system (1). 

Consider a wave equation of the type (29) 

#-il>xx + W{x)i> = 0, (37) 

with an arbitrary potential W(x) (whose specific example is the above potential V c q). If there is a 
function S(x) such that W(x) is presented in the form 

W(x) = S 2 (x) + S x , (38) 
then Eq. (37) is rewritten as follows: 

^ + (d x + S)(-d x + S)^ = 0. (39) 
Now, if we introduce the new function 

X = (~d x + Sty, (40) 

then, applying the operator — d x + S to the left-hand side of Eq. (39), we obtain the following wave 
equation for X'- 

X~Xxx + W rcg (x)x = 0, (41) 
with the new effective potential 

W reg (x) = -S x + S 2 = -W(x) + 2S 2 . (42) 
If a static solution ip s (x) of Eq. (37) is known, so that tp SjXX = W(x)ip s , then we can choose 

S(x) = ^ s , x /i> s (43) 

to carry out the above transformation. 

Generically, the function U in (24) and hence the potential (30) behave near a throat as r'~ 2 ~ 
(u — u th )~ 2 ~ x~ 2 , where, without loss of generality, we put x = at the throat. Assuming that 
the potential W(x) behaves in such a way, let us look if a transition to W rcg can really remove this 
singularity. Above all, we see that according to (42), such removal is only possible if W — > +oo as 
x — > since we are dealing with real quantities. Thus a potential wall in W(x) can be removed 
but a potential well cannot. 

A positive pole W ~ x~ 2 can be removed in W reg if 

S »j -l/x =^ oc 1/x. (44) 

It is a necessary condition for regularizing the potential. Besides, to avoid a singularity of Wreg, 
ip s must be nonzero in the whole range of u. Moreover, according to (42), it is clear that near the 
throat x = in this case 



W(x) w 2/x 2 . 



(45) 
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Thus, to be regularized by the procedure described, the pole in W(x) must behave as (45). Let us 
show that this is generically the case for the potential W = V e s(x) given by Eq. (30). 
Suppose such a generic situation, so that 

(i) the function f3(x) in the background metric is expanded near its minimum (the throat) in powers 
of x as follows: 

/3( x ) = f3 + \f3 2 x 2 + |/3 3 x 3 + . . . . (46) 
where A),2,3 are constants; 

(ii) the background quantity (f>'(x) ^ at x — 0. 

Here and till the end of the section, we use the coordinate freedom to choose the "tortoise" 
radial coordinate u = x specified by the condition a — 7. All functions are considered as power 
series in x at small x. 

Let us estimate V e s (30). The term that determines the pole at x = is 

W po Ux)=e^(V-e~^, (47) 

where we have put e = — 1 since throats are possible only in this case. By (46), 

& = fax + \fox 2 + .... (48) 

Now we use the equations governing the static configuration: from (15) it follows that e 27 (e _2/3 — 
V) = (3 XX at x = 0, and then from (16) we find that $\ = 2f3 xx at x — 0. Thus from assumption 
(ii) it follows /3 2 7^ 0. Substituting all this to (47), we find that it behaves precisely as required in 
(45). 

Thus we have shown that, for a (generic) throat in a solution to Eqs. (12)— (16), the effective 
potential V e & for spherically symmetric perturbations satisfies the necessary condition for regular- 
ization by the above method. 

Whether or not this regularization really works and leads to a regular boundary-value problem 
for the perturbations, should be investigated for specific background configurations. (In particular, 
one should also take into account the other singular term in the potential V e s, proportional to 
V(j>/ fix'- the terms oc 1/x evidently depend on the finite part of S(x).) 

A positive example of such a study, concerning anti-Fisher wormholes, is known from [18]; we 
here show that the other two branches of the anti-Fisher solution are also such examples. 



2.6. Regular perturbations near a throat 

Suppose we have found a solution x{ x i t) to Eq. (41), satisfying the appropriate boundary conditions. 
The function x is regular at x = since the potential W rcg is regular there. If x is a growing 
function of t, it probably indicates an instability of the initial static configuration; but it is indeed 
the case only if this x( x ,t) creates regular perturbations of the metric functions a, (3, 7 and the 
scalar field 0. 

Let us look how it happens. Given x( x , t) , a solution to (37), or (29), is found as ip — (d x + S)x- 
Generically, x is finite at x = while S ~ 1/x at small x, hence ip ~ 1/x, and according to (28) we 
obtain 5<p — > 00 at the throat. This result is in fact quite natural since the relation (28) corresponds 
to the gauge 5r = 0, in which the throat radius is fixed, whereas we were seeking perturbations 
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with nonzero 5r on the throat. So it is necessary to pass on to another gauge, which is easily done 
due to gauge invariance of the quantity ip given by 

i/j = r5<j> - —5r. (49) 

Namely, a finite expression for 5r is obtained in the gauge 5(f) = provided <f> x (0) ^ since then 

5r = -— t/j, (50) 
r(p x 

while the product r x ip is finite. It remains to find 5a and £7 from the perturbation equations in 
the gauge 5(f) = 0. From Eqs. (18) and (17) we find 

(5 x 5a = 5f3 x + 5(5((5 x - lx ), (51) 
5l x = 5a x - 25 /3 X - ^ e 2a 5a, (52) 

<Px 

but here we are again facing a problem: according to (51), in general 5a diverges at the throat where 
f3 x = 0. This divergence is only avoided if the right-hand side of (51) behaves like (3 X ~ u — u t h ~ x. 

Surprisingly, it is the case in a generic situation, as can be verified in a general form using 
near-throat expansions. Indeed, let us preserve the above assumptions (i) and (ii) and assume, in 
addition, that 

(iii) the function x( x jt) that solves Eq. (41) is finite and nonzero at x — 0, 

(iv) the function S(x) behaves at small x according to (44). 

Our task is to estimate the right-hand side of Eq. (51) in the order 0(x°) : if it is zero, it means 
that (5 x 5a ~ x ~ (3 X and thus 5a (0) is finite; the remaining metric perturbation £7 is then found 
from (52) and is also finite. 

Taking 5(3 = 5r/r from (50) and substituting ip as ip = Xx + Sx ~ ~Xo/ x , where xo — x(0)> 
we find 



Then we substitute this expression to (51) to obtain 

P x 5a » [ - <p xx (3 2 + \<j> x {p 3 - 2 lx (3 2 ) 

X L 

where all quantities are taken at x — . Now, <p xx can be expressed from the background equation 
(12), (3 XX from (15); we can use the fact that f3 2 = /3 XX (0) etc., and we can also ignore all terms 
proportional to (3 X . After these substitutions we finally obtain 



(3 x 5a oc e^Vs 



+ ¥x 



l j2 

. X =0 



But the expression in the square brackets vanishes due to the difference of Eqs. (13) and (14), which 
proves that (3 x 5a = 0(x) and thus 5a(0) is finite. 

We conclude that under the generic assumptions (i)-(iv), regularization of the potential V e s 
always leads to finite perturbations of the background static solution. 
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3. Instabilities of the Fisher and anti-Fisher solutions 
3.1. The static solutions 

Let us recall the well-known static, spherically symmetric solutions to the field equations (2) and (3) 
for zero potential, V = 0. In the case e = +1 this solutions was found by I.Z. Fisher in 1948 [8] and 
afterwards repeatedly re-discovered. For e — — 1 the corresponding solution was first obtained, to 
our knowledge, by Bergmann and Leipnik [28] and also repeatedly re-discovered. But these authors 
used the curvature coordinates [i.e., the condition u = r in terms of the metric (4)], which are not 
well suited for the problem, and maybe therefore they did not give a clear interpretation of the 
solutions. 

The solution can be written jointly for e = ±1 if one uses the harmonic coordinate u = v in 
the metric (4), corresponding to the coordinate condition a(v) = 7(f) + 2/3(v) [10]. The scheme 
of obtaining it is as follows. In these coordinates, the scalar field equation and the Q component 
of (3) read simply 0" = and 7" = 0, respectively (the prime here stands for d/dv), so that, 
choosing the zero point of and the time scale, without loss of generality we can write = Cv 
and 7 = —mv with C,m — const . Furthermore, a sum of the components (°) and Q) of (3) leads 
to the Liouville equation (j3 + 7)" = e 2/3+27 whose solution is 

{k" 1 sinh kv, k > 
v, k = (53) 

k" 1 sin kv, k < 0. 

As a result, the solution as a whole reads 



2mv r J 2 

d S 2 =e~ 2mv dt 2 ' 



s 2 (k, v) ls 2 (k, v) 



+ dn 2 



Cv, (54) 



where the integration constants m (the Schwarzschild mass), C (the scalar charge) and k are 
related by the equality 

2k 2 signk = 2m 2 + eC 2 , (55) 

obtained after substituting all found functions to the Einstein equation G\ — . . . which does not 
contain second-order derivatives in v and is an integral of other equations. 

The coordinate v is defined in the whole range v > for k > and in the range < v < ir/\k\ 
for k < 0. The value v = in all cases corresponds to flat spatial infinity, so that at small v the 
spherical radius is r(v) ~ 1/v, and the metric becomes approximately Schwarzschild. 

In the case k > , it is helpful to pass over to the quasiglobal coordinate coordinate u (defined 
by the condition a + 7 = in (4)) by the transformation 

e -2kv = 1 _2k/ u =: P(u), (56) 
and the solution takes the form 

ds 2 = P a dt 2 - P~ a du 2 - P l - a u 2 dtf, = ~ In P{u), (57) 

ZK 

with the constants related by 

a = m/k, a 2 = l-eC 2 /{2k 2 ). (58) 

The Fisher solution [8] corresponds to e — +1, hence according to (55), it consists of a single 
branch k > and, in (57), \a\ < 1. It is defined in the range u > 2k, and u = 2k is a naked central 
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(r = 0) singularity which is attractive for m > and repulsive for m < 0. The Schwarzschild 
solution is restored at C = 0, a = 1 for m > and at C = 0, a — — 1 for m < 0. 

The solution for e = — 1 (that is, for a phantom scalar field) is conveniently termed the anti- 
Fisher solution, by analogy with de Sitter and anti-de Sitter. According to three variants of the 
function (53), this solution splits into three branches with the following properties. 

Branch A, k > 0: the solution again has the form (57), but now \a\ > 1. For m < 0, that is, 
a < —1, we have, just as in the Fisher solution, a repulsive central singularity at u = 2k. The 
situation is, however, drastically different for a > 1 . 

Indeed, the spherical radius r in all such cases has a finite minimum at u — u t h = (a + l)k, 
corresponding to a throat of the size 

r(«th) = r th = k(a + l)( a+1 )/ 2 (a - (59) 

and tends to infinity as u — > 2k. Moreover, for a = 2, 3, ... the metric exhibits a horizon of order 
a at u = 2k and admits a continuation to smaller u [12]. A peculiarity of such horizons is their 
infinite area. Such asymptotically flat configurations with horizons of infinite area have been termed 
cold black holes (CBHs) [11] since all of them have zero Hawking temperature. The throat radius 
(59) does not coincide with the Schwarzschild mass m = ak (in usual units, half the Schwarzschild 
radius 2Gm/c 2 ) but is of the same order of magnitude. 

Furthermore, it can be verified using (6) [12] that the metric (57) has a curvature singularity 
at u = 2k if a < 2 (except for the Schwarzschild case a = 1), a finite curvature if a = 2 and zero 
curvature if a > 2. 

For non-integer a > 2, the qualitative behavior of the metric as r — > 2k is the same as near 
a horizon of infinite area, but a continuation beyond it is impossible due to non-analyticity of the 
function P a (u) at u — 2k. Since geodesies terminate there at a finite value of the afline parameter, 
this is a space-time singularity (a singular horizon as it is named in [11]) even though the curvature 
invariants tend there to zero. 

Branch B, k — 0: the solution is defined in the range u E K + and is rewritten in terms of the 
quasiglobal coordinate u — 1/v as follows: 

ds 2 = e~ 2m/u dt 2 - e 2m/u [du 2 + u 2 dVl 2 ], <f) = C/u. (60) 

As before, u = oo is a flat infinity, while at the other extreme, u — > 0, the behavior is different 
for positive and negative mass. Thus, for m < 0, u = is a singular center (r = 0), while for 
m > 0, r — > oo and all Ki — > as u — > 0. This is again a singular horizon: despite the vanishing 
curvature, the non-analyticity of the metric in terms of u makes its continuation impossible. The 
throat occurs at u — m and has the size e • m , e being the base of natural logarithms. 

Branch C, k < 0: the solution describes a wormhole with two flat asymptotics at v — and 
v — n/\k\. The metric has the form [9, 10] 



ds 2 = e~ 2mv dt 2 k 



2 D 2mi) 



k 2 du 2 



+ dQ 2 j 



sin (kv) Lsin (kv) 

= e~ 2mv dt 2 - e 2mv [du 2 + (k 2 + u 2 )dn% (61) 

where v is expressed in terms of the quasiglobal coordinate u, defined on the whole real axis, by 
kv = cot^ 1 (u/k) , where we have denoted — k = k > 0. Ifm>0, the wormhole is attractive for 
ambient test matter at the first asymptotic (u — > oo) and repulsive at the second one (u — > — oo), 
and vice versa in case m < 0. For m = one obtains the simplest possible wormhole solution, 
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sometimes called the Ellis wormhole, although Ellis [9] actually discussed these solutions with any 
m. 

The wormhole throat occurs at u = m and has the size 

r th = ( m 2 + P)V2 exp ( ™ cot-i ^ . ( 62 ) 

\k k J 

3.2. Perturbations: the Fisher solution 

For a massless scalar field [V = in (1)], the effective potential in (29) takes the following form in 
terms of the quasiglobal coordinate u: 

A(f)' 2 Ar" A'r' A<p' 2 A A 2 r' 2 , x 

V « = + — + — = " e ^T + -2 ~ (63) 

where A(u) := e 27 ^ = e~ 2a ( u \ and the second equality in (63) follows from Eq. (15). 

Calculating V e R for the solution (57), we find a common expression for both e = +1 and e = — 1 : 

~ j. / \ , n2a 2cm 3 - 3(1 + a) 2 ku 2 + 2(3 + 4a + 3a 2 + 2a 3 )k 2 u-(l + a) 4 k 3 

V ^ = kP u 2 (u-k(l + a)) 2 (u-2k) 2 • (64) 

Since in the Fisher solution a < 1 , the binomial u — k(l + a) is positive at all u > 2k , and the only 
singularity in V c g is u — > 2k, coinciding with the singularity of the background solution. Near the 
singularity, at which according to (27) we can put x — 0, V c q{u) ~ — l/(4x 2 ), a negative pole in 
agreement with [26] and many subsequent papers. 

The boundary condition at spatial infinity (u — > oo, x — > oo) is natural: 5(f) — > 0, or ip — > 0. 
For u — > 2k, where the background field tends to infinity, the boundary condition is not so 
evident. In [26] and other papers, dealing with minimally coupled or dilatonic scalar fields, the 
minimal requirement was used providing the validity of the perturbation scheme: 

\5<p/<f)\ < oo. (65) 

(The requirement of absence of ingoing waves then does not lead to further restrictions.) Under this 
boundary condition it is easy to conclude that there are solutions to the Schrodinger-like equation 
(32) with any u 2 < 0, which means that the static field configuration is unstable, in agreement 
with the previous work [26] (see also [12] for details). 

3.3. Perturbations: the anti-Fisher solution 
Branch A 

The effective potential has the same form (64), but now, since a > 1 (we restrict ourselves to this 
case providing m > 0), the potential has a positive pole at u — u t h = k(a + 1) > 2k, the throat in 
the background configuration. 

To remove this singularity, we were able to find the following simple static solution to Eq. (29): 

ip s +{u) oc r(u)— — — , (66) 
u - u th 

Applying the technique described in Sec. 2.5 with these ip s (u) , we the new effective potential W rcg = 
W A {u): 



W A (u) =1 



2k\ 2a N(u) 



u J u 2 {u — 2k) 2 (au — u th ) 2 ' 
N(u) = 3(1 + a) 4 k A - 2a(6 + 19a + 16a 2 + 3a 3 )A; 3 M 

+ 3a 2 (9 + 10a + a 2 )A;V - 2a 2 (4 + 5a)A;w 3 + 2aV. (67) 
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Figure 1: The regularized effective potential Wa (left panel) and the time-domain profiles (right 
panel) for the Branch A solutions with a = 3/2 (blue), a = 2 (green), a = 3 (red). Smaller values 
of a correspond to deeper potential wells and a more rapid growth of perturbations. 




Figure 2: The potential Wa (left panel) and the time-domain profile (right panel) for the Branch 
A solution with a = 1 . 



The potential (67) has no singularities at r > 2k and is partly negative, which should in 
general lead to an instability. To prove it, one can use the method of time-domain integration [29] 
allowing for following the time evolution of the perturbations under prescribed initial and boundary 
conditions. The latter, as usual, must provide the absence of ingoing waves from the boundary, 
and in our case it is sufficient to simply require ip — >■ and x ~ ?• as x — )• ±oo for all three 
branches of the anti-Fisher solution. The reason is that the effective potential V e g as well as the 
regularized potentials vanish at large \x\ , and therefore modes of interest, those with u 2 < 0, should 
exponentially decay at large 

Examples of plots for the potential (67) and the results of time-domain integration are shown in 
Fig. 1. By fitting of the profile we find that the perturbations grow approximately as ip oc e a25 */ m . 

In the limit a — > 1 the regularized potential still has a negative gap, however, at a = 1 no 
growing mode is observed, and a stationary solution dominates at late times (Fig. 2), while at 
any finite a — 1 the perturbation does grow. This shows how the instability is "dying out" when 
approaching the Schwarzschild solution. One can recall that in the genuine Schwarzschild case there 
is no scalar field, and the modes we are considering here simply do not exist. 
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Figure 3: The potential Wb (left panel) and the time-domain profile (right panel) for the Branch 
B solution. 



Branch B 

For the solution (60), the effective potential (63) 

/ 4m \ 2m 3 — 3mu 2 + 4m 2 « — m 3 
Vtfr(u) = mexp —-^ 



is singular at u t h = m. We again find a static solution to (29) 



u 2 e m/u 



u — m 



and perform the transformation described in Sec. 2.5. 

The regularized effective potential has the form of an inverse potential barrier: 

, . / 4m \ 3m 2 — 10mw + 2u 2 
W B {u) = exp , 

\ U J u 

The perturbations grow approximately as ip oc e°' 23t//m (see an example in Fig. 3). 
Branch C 

For the wormhole solution (61), the potential (63) has the form 



V eS (u) 



_ imv 2k A + k 2 (3m 2 — 2mu + 3m 2 ) — m(m 3 — Am 2 u + 3mu 2 — 2u 3 ) 

(m — u) 2 (k 2 + u 2 ) 2 



(6J 



(69) 



(70) 



(71) 



It again exhibits a positive pole on the throat. To regularize it, we can take any of the two static 
solutions 



/ / \ T'(li] r 2 o /t 2 \i 

ip s —\\u) oc [k +m — mv{k +mu)\, 



u — m 



, x r ( U ) /,2 \ 

ip s -2{u) oc [k + mu) 



u — m 



(72) 



(73) 



The first of them coincides with the solution found by Gonzalez et al. in [18]. As a result, we 
reproduce their regular potential 



Wcx{u) = e 



-4m» 



NAu) 



(k 2 + u 2 ) 2 (k 2 (mv — 1) + m 2 (uv — l)) 2 ' 



(74) 
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Figure 5: The potential Wei (left panel) and the time-domain profile (right panel) for Branch C 
(k = m = 1). 



where 

Ni(u) = -3k 6 (mv - l) 2 + m 4 (3m 2 - lOmu + 2u 2 )(uv - l) 2 

+ k 2 rr?{uv — l)(8u + 2m 2 v + m(5 — 15uv)) + k 4 m(mv — l)(2u + m 2 v + m(9 — 12m')). 

A further investigation revealing an unstable mode of perturbations is described in detail in [18] 
(Fig. 4). 

The second static solution (73) can be used to demonstrate that, despite the different regularized 
potentials, the physical result, namely, the perturbation growth rate remains the same. Indeed, the 
regularized potential produces by the solution (73) reads 

rrr / x Am, ,~ 3 & 6 + k^m (tti — 1 2u) + ^tti 2 (2m — 15m) u + m 2 u 2 (3m 2 — 10mw + 2u 2 ) ) , . 

W C2 (u) = e imv ^ '- j—± ^ ^, 75 

[k z + muyik* + u A y 

and leads to perturbations growing as shown in Fig. 5, at the same rate as in Fig. (4). 



4. Concluding remarks 

We have demonstrated the instability of all anti-Fisher solutions to the Einstein-scalar equations 
under spherically symmetric perturbations, thus confirming and extending the conclusions of [18] 
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made for Branch C solutions (wormholes) . It turns out that in almost all cases the characteristic 
time of perturbation growth is of the order of the time needed for a light signal to cover a distance 
equal to the throat radius. 

We have found out that the S-deformation method of regularizing the effective potential for the 
perturbations, having a positive pole at a throat (a minimum of the spherical radius r(u)), used 
in [18], is applicable to any static, spherically symmetric configurations of self-gravitating scalar 
fields with self- interaction potentials. We have shown in a general form that, under some generic 
assumptions [items (i)-(iv) in Sections 2.5 and 2.6], the potential V e g near the throat has the form 
admitting regularization, and that a regular mode found as a solution to the regularized equation 
leads to a regular perturbation of the initial configuration. The latter circumstance is quite non- 
trivial because regularity is required not only for solutions to the master equation but also for all 
metric coefficients. 

However, application of this methodology to specific solutions with self-interaction potentials 
V{4>) 7^ faces the problem of explicitly finding the function S(x) satisfying Eq. (38), or, equiva- 
lently, a proper static solution to the master equation (a zero mode). We hope to extend this study 
to some configurations with nonzero V((f>) in the near future. 

One more important extension of the presently obtained results should be mentioned: Eqs. (2), 
(3) are well known to be not only the equations of general relativity but also the Einstein-frame 
equations of the general scalar-tensor and f(R) theories of gravity, see, e.g., [30-33]. Note that in 
f(R) theories, the Einstein frame always contains a scalar field with nonzero V((f>). The observable 
physics in such theories is usually described in the corresponding Jordan conformal frame, whose 
metric differs from that of the Einstein frame only by a conformal factor that varies from theory 
to theory. A transition from one frame to the other is, from the viewpoint of differential equations, 
simply a substitution, therefore the stability study can safely be performed in the Einstein frame. 
Then, if the conformal factor is everywhere regular and nonzero (including the end points of the 
coordinate range), the conformal mapping preserves the boundary conditions for regular pertur- 
bations, and the stability conclusions obtained in the framework of general relativity are readily 
extended to the corresponding solutions of these generalized theories. In particular, we can assert 
that all vacuum static, spherically symmetric solutions in scalar-tensor theories, connected with the 
(anti-)Fisher solutions by everywhere regular conformal factors, are unstable. 

However, in many cases the conformal factors bear nontrivial features, i.e., somewhere blow up 
or vanish, and this can affect the boundary conditions for the perturbation equations; in any such 
case a separate study is necessary. Consider, for instance, the counterpart of Fisher's solution (57) 
in the Brans-Dicke (BD) scalar-tensor theory, where the Jordan-frame metric is g 3 = (l/$)p MI/ ; 
$ = exp(0/ s/uj + 3/2) is the BD scalar field, oj > —3/2 is the BD coupling constant, and |a| < 1. 
It is the so-called Brans class 1 solution. In both frames, the value u = 2k is a naked singularity 
(see more details on these solutions in [11,34]). In the stability study, to formulate a boundary 
condition at this singularity, for Fisher's solution we have used the minimal requirement (65), 
\8(j)/<j)\ < oo, providing the validity of the perturbation scheme, and we then concluded that the 
background solution is unstable. In the BD picture, it is more reasonable to require that the 
perturbed conformal factor l/<3> behave not worse than the unperturbed one, i.e., |5$/$| < oo. 
However, since 5$/$ ~ 5(f), we arrive at the condition \5<f)\ < oo which is more restrictive then (65) 
if — > oo, and this made us conclude in [11] that this BD solution is stable. 

Another example of using the conformal mapping between Jordan and Einstein frames for sta- 
bility studies can be found in [22,23], where the instability was proved for electrically neutral 
and charged wormholes supported by nonminimally (in particular, conformally) coupled scalar 
fields [10,23,35]. In this case, there is a drastic difference between the manifold structures in the 
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two conformal frames: in the Einstein frame, without an electric charge, it is the Fisher solution 
with the metric is that has a singularity at u — 2k. In Jordan's, this singularity is removed due 
to the conformal factor, the solution is continued beyond this (now regular) sphere and has a flat 
asymptotic at the other end. Such wormholes proved to be unstable [22,23], but their instability is 
of quite different nature than that described in this paper and in [18]: it is related to a negative pole 
of the effective potential (64) at u = 2k for Fisher's solution (57), e = +1 and a similar singularity 
in its counterpart with an electric charge. 

These examples show that there remains quite a lot of work in studying the stability of more 
complicated solutions of scalar-tensor and f(R) theories. 
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